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Abstract. - To enlighten the nature of the order-disorder and order-order transitions in block 
copolymer melts confined in narrow capillaries we analyze peculiarities of the conventional Landau 
weak crystallization theory of systems confined to cylindrical geometry. This phenomenological 
approach provides a quantitative classification of the cylindrical ordered morphologies by expan- 
sion of the order parameter spatial distribution into the eigenf unctions of the Laplace operator. 
The symmetry of the resulting ordered morphologies is shown to strongly depend both on the 
boundary conditions (wall preference) and dimensionless parameter q*R, where R is the cylinder 
radius and is the wave number of the critical order parameter fluctuations, which determine 
the bulk ordering of the system under consideration. In particular, occurrence of the helical mor- 
phologies is a rather general consequence of the imposed cylindrical symmetry for narrow enough 
capillaries. We discuss also the ODT and 00T involving some other simplest morphologies. The 
presented results are relevant also to other ordering systems as charge-density waves appearing 
under addition of an ionic solute to a solvent in its critical region, weakly charged polyelectrolyte 
solutions in poor solvent, microemulsions etc. 



Introduction. — During the last years there has been 
large interest in the behavior of block copolymer (basically 
diblock copolymer) melts confined in narrow capillaries, 
which were studied both experimentally [1-6] and via com- 
puter simulation [7-9] as well as by numerical calculations 
within the self-consistent field theory [1,10,11]. The most 
characteristic of these morphologies are lamellae ordered 
normal to the cylinder axis z (slab morphology [10] ) and 
coaxial shells or multiwall morphologies [7, 10] as well as 
various helices [1,8,9], the morphologies' symmetry being 
strongly depended on the ratio of the gyration radius Rq 
of the diblock macromolecules under consideration and the 
radius R of the cylinder the macromolecules are confined 
to. However, some important issues have not been ad- 
dressed yet. In particular, the various morphologies (up 
to 20 in ref [11] and up to 30 in ref [9]) are distinguished 
only visually by their snapshots and no attempts to relate 
them to a quantitative order parameter are made until 
now. Meanwhile, without a quantitative definition of an 
order parameter and, accordingly, finding the quantitative 
change of such an order parameter at the phase transitions 



one cannot discuss the transitions between the morpholo- 
gies rigorously. Probably, this is why no accurate phase di- 
agram of the order-disorder and order-order transitions for 
cylindrically confined block copolymer melts is presented 
until now. In particular, the phase diagram presented in 
ref [11] has two important shortcomings: i) it is based on 
a visual distinguishing of the morphologies only, and ii) 
the authors assume no modulation of the block copolymer 
composition along the z-axes, which, obviously, excludes 
the slab and helical morphologies. Besides, an alternation 
of some coexisting phases is to occur in cylindrical capil- 
laries already by virtue of the Landau theorem [12] as in 
any effective ID system. 

Another shortcoming of all the aforementioned stud- 
ies [1-11] is that they are specifically block copolymer 
oriented whereas ordering in cylindrically confined sys- 
tems seems to be much more general. Indeed, the very 
existence of the various morphologies in block copoly- 
mer melts is due to the instability of their homogeneous 
phase with respect to spatial fluctuations of the polymer 
concentration with a finite spatial period [13-15]. Thus, 
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the ODT in block copolymers is only a particular case 
of the well known weak crystallization instability [16, 17]. 
The Landau weak crystallization theory (WCT) provides 
a common physical background for describing such seem- 
ingly different physical phenomena as i) forming weakly 
segregated morphologies in molten block copolymers [14], 
therewith the polymer specific features appear at the stage 
of the microscopic calculation of the Landau expansion co- 
efficients only, ii) appearance of the blue phases in liquid 
crystals [18,19], Hi) charge-density waves generation upon 
addition of an ionic solute to a solvent in its critical region 
[20,21], iv) microphasc separation in weakly charged poly- 
electrolyte solutions [22,23] and v) microemulsions [24]. 

So, the peculiarities found for cylindrically confined 
block copolymer melts are expected to exist for other 
nano-ordered systems either and the aim of the present 
Letter is to present a phcnomcnological theory of the or- 
dering in narrow capillaries based on the general weak 
crystallization paradigm. 

1. We start with the standard expression for the free 
energy of the weakly crystallized systems described by a 
non- uniform scalar order parameter profile <&(]"), which 
without any loss of generality could be written as the Lan- 
dau expansion in powers of <f> up to the 4th order: 

V -^= [ (r$ 2 (r) + ((l + ^ 2 A)$(r)) 2 (1) 



+a$ 3 (r) + $ 4 (r)) dr, 

where v is a constant having dimensionality [Z 3 ] to be de- 
fined from microscopic considerations and r is an effective 
dimcnsionless temperature measured from the instability 
boundary for bulk. In bulk the quadratic term in (1) reads 



F (2) 
bulk 



r + (i-<z 2 A* 2 ) 2 ' 



(2) 



where $ q = J $ (r) exp (iqr) d r is the Fourier component 
of the order parameter. Obviously, for r > a mini- 
mum (at least, a metastable one) of the free energy (1) 
is provided for the disordered state ($(r) = 0), whereas 
for r < the disordered state is absolutely unstable 
with respect to the growth of the Fourier components $ q 
with wave numbers q — |q| close to the value q = q*. 
Thus, the parameter q% > characterizes the wave length 
L = 2tt/ g* of the critical order parameter fluctuations, 
which destroy the disordered state in bulk. 

The coefficient a appearing in the free energy expression 
(1) depends on the physical nature of the system. In this 
paper we restrict ourselves to the case a = 0, when the 
only ordered phase appearing in bulk systems described by 
(1) is lamellar, and show that even in this case the presence 
of confinement could result in various morphologies formed 
by straight and curved rod-like units. 

2. In the presence of a confinement the order parame- 
ter ^(r) is to be expanded not in the Fourier harmonics, 
which are the cigenfunctions of the Laplace operator in the 



infinite volume, but in the eigenfunctions of the Laplace 
operator 

A* (r) = -EV (r) (3) 

for the specified confinement geometry. In the cylindrical 
co-ordinates, which are the natural choice for ordering in 
capillaries, the eigenvalues and the corresponding eigen- 
functions read 

oo oo 



E = a 2 m , n +p 2 , tt(r)= E 



m— — oo n— 1 



f 



dp 



A P! m,n exp (ipz) exp (irrup) m ,„ (r) — (4) 

Zir 



Here the order parameter modulations along the cylin- 
der axis and in polar angle are characterized by the wave 
number p and the integer number m, respectively, whereas 
<t>m,n{ r ) are the eigenfunctions of the 2-dimensional radial 
Laplace operator: 



d4 

r dr Y dr 



m 



(5) 



4>m,n{r) are the Bessel functions [25] J m (a m ,n f)i the 
eigenvalues a 2 nn being determined by the corresponding 
boundary condition. Choosing the reflecting boundary 
condition 

d*/dr\ r=R = 0, (6) 

which is often used to describe the behavior of confined 
incompressible block copolymer melts [26, 27], we get 
a m n = s m nl ^ i where s TOi „ are the locations of extrema 
of the Bessel function J. m (x) numerated in ascending or- 
der. 

3. Substituting eqs (4) into (2) we get 



/oo 00 00 
2 A m,n(T,P,p)Mn,n \Ap,r 
-oo m=-oo n=l 



dp, 



A m ,„(r,p,p) = r + (p 2 /q 2 - n 2 m . n ) 2 , (7) 
where p = q*R is the reduced capillary radius, W m , n = 
Jo Jm { s m,n %) x dx and 

K 2 m ,n = 1 - S m , n /P 2 - (8) 

The condition providing instability of the disordered 
state with respect to infinitesimal fluctuations of the 
(m, n)-mode is 



A m p)<0, 



(9) 



where p* m n is location of the absolute minimum of the 
function (7). 

It follows from eqs (7), (8) that the condition (9) takes 
two rather different forms: 

P*m.n = q* K rn,n, T < 0, for S m . n < p (10) 

P*m.n = 0, T < -«4 )Tl , for S m , n > p. (11) 
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Fig. 1: A typical succession of the ordered (black) and dis- 
ordered (white) capillary segments for a finite value of the 
Ginzburg parameter. 



In other words, given a finite reduced capillary ra- 
dius p all unstable modes with s m ^ n > p result in an 
angular and radial modulation only and their contribu- 
tion to the total order parameter profile is proportional 
to Jm{s m ,nT / R) cos(mtp). On the contrary, the unstable 
modes with s myn < p possess both an axial and, generally, 
angular as well as radial modulation. 

In the narrow capillaries obeying the inequality 



p < si,i = 1.841 



(12) 



the only mode (0, 1), which corresponds to the cigcnhmc- 
tion 

$o,i (r) = A cos(q*z) J (s 01 r/R) = A Q cos(q*z), (13) 

is axially modulated (s ,i = 0). Within the temperature 
interval (— < t < 0) minimizing the free energy (1) 
we get 



Al = -2t/3, F bulk = -T{L/R) r 2 /Gi, 



(14) 



L and Gi = 6v/(7ri? 3 ) being the length of the capillary 
and Ginzburg parameter, respectively. 

4. Indeed, by virtue of the Landau theorem [12], the or- 
dered phase appears as an alternating succession of the 
ordered and disordered segments of various, generally, 
lengths (see fig. 1) rather than ordering in the whole capil- 
lary. Accordingly, the capillary partition function could be 
written similarly to that of Id Ising (helix-coil transition) 
model [28,29]: 



Z(L) 



exp (pL) dp, 



(15) 



where the Laplace transform Z{p) reads 



Z(P) = f 
Jo 



exp (-pL) Z(L) dL = aiajZijip), 



\Zi 3 (p)\\ = 



p —a 
-a p + f 



(16) 



where Tf is the specific (per unit length) free energy of the 
ordered phase, i,j = 1,2 and a, at, oi are, respectively, 
the weights to be assigned to transient zones between the 
ordered and disordered segments as well as those between 
the end-walls of the capillary and the disordered and or- 
dered segments. 1 It follows from eqs (15), (16) that the 



1 Strictly speaking, the free energy F(l) of a capillary segment 
of the length / is slightly different from fl since for finite segments 
one should use Fourier series and ^-dependent boundary conditions 
on the end- walls rather than Fourier integrals (5b), which difference 
is neglected in this paper. 



length fraction of the ordered segments is a continuous 
function of the parameter u = f /2a (u < for r < 0) 



^ord = -L- X d In Z(L)/df= (Vu 2 + l-u)/( 2^ + 1) 

(17) 

Taking into account that a ~ R 1 exp(fR) and, therefore, 
|u| <~ \fR\ exp \ fR\, we conclude that capillary ordering is 
not a genuine phase transition but a co-operative one, the 
location and width of the transition region being r = 
and \t\ — (Gi) 1 / 2 . Thus, if the inequality (12) holds and 
Gi <C 1 then, as is seen from eq (14), for practical purposes 
the line t = could be considered as the 2nd order ODT 
line. 

The presented consideration of capillary ordering is ex- 
tended straightforwardly to the case when the free energy 
(1) has N > 1 minima corresponding to different mor- 
phologies, which provides the set /i,../jv of the specific 
(per unit length) free energies. Accordingly, the matrix 
Zij and vector Oi appearing in eq (16) become N + 1- 
dimcnsional ones. Then an equilibrium between the seg- 
ments with different ordered morphologies and the dis- 
ordered segments would occur, character of which would 
change on crossover lines. In the limit Gi — > these 
crossover lines transform into the 1st order transition lines 

mm(f 1 ,..f N ) = f A = f B (18) 

between the ordered morphologies A and B. 

5. One more mode, (1,1), is axially modulated for 

si.i < p < s 2 ,i = 3.054, (19) 

two types of the axially modulated profiles being possible: 



(20a) 



(r) = A h cos {<p ± Pl>1 z) Ji (s M r/R) , 
$g (r) = A s cos {ip) cos {p 1A z) J x (s M r/R) . (20b) 

Further we refer to the profiles (20a) and (20b) as the run- 
ning along the axes z (helical) and standing waves, respec- 
tively. The specific free energies obtained by substituting 
the order parameter profiles (20) into the free energy (1) 
read: 



fh = TM m ,nA 2 h 



A 4 

U = (r/2)AA m ,„A 2 + (9/32)S 



(3/4)B m ,„^, 



(21a) 
(21b) 



where / = fRG± and B m . n = J\{s m , n x) x dx, (the 
difference in the corresponding numerical coefficients is 
due to difference in averaging the powers of cosines (for 
Fh) and the products of such powers (for F s )). Minimizing 
the free energies (21) with respect to the corresponding 
amplitudes results finally in the expressions 

f h = -2T 2 J>C, n B-] n , / s = -(4/3)rX,„C. (22) 

As is seen from (22), the single helical wave is always more 
thermodynamically advantageous than the standing one 
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just due to the symmetry properties. The free energies of 
the left- and right-hand helical waves are identical and the 
capillary would split into the left- and right-hand helical 
segments with the average segment length ~ er~ 2 if the 
helical state is dominant. 

For p > s 2 ,i at least one more mode (2,1) becomes axi- 
ally modulated, the contribution of the mode into the total 
order parameter profile resembles the double (DNA-like) 
helix: 

(r) = A 2h cos (2ip ± P2.iz) J 2 (82,i r/R) (23) 

To check whether and when one non-uniform mode gives 
way to another more thermodynamically advantageous 
one or a mixture of two modes with further change of 
the temperature and capillary radius, we substitute into 
the total free energy (1) a 2- mode trial order parameter 
profile 

<E>(r) = AiM cos(g*K mi , ni z - mi<j>) J mi (s mi , ni r/R) (24) 
+A 2 N 2 cos(q*K m2 ^z - m 2 <j>) J m2 (s m2 ,„ 2 r/R) 

The proper choice of the normalizations A/i leads to the 
reduced specific free energy 

7=3minF(A 1 ,A 2 ), F(A 1 ,A 2 )=t 1 A 2 +t 2 A 2 
+ (3/8)(C lCll 4 + 2c 12 A 2 A 2 + C 2 c 22 A 4 2 ), (25) 

where = 2 J* J^. (s mi , ni x) J^. (s mj , nj x) xdx j A/iA/}, 

cu = Jo J„. (s muni x) xdx j N 2 , n = t + «4 i)fli for the 

modes (10) and r» = r for the modes (11), Cj = 2/3 if the 
i-th mode is a mode (0, n) belonging to the type (10) and 
d = 1 otherwise. 

The straightforward analysis shows that the minimum 
of the free energy (25) is achieved in a 1-mode state only 
(Ai = or A 2 = 0) if the determinant | is negative, the 
transition between the modes being the discrete 1st order 
phase transition in the limit Gi — > 0. If \cij\ > then 
the minimum of the free energy (25) could correspond to 
a 2-mode state (A\ ^ and A 2 ^ 0), the second mode ap- 
pearing (in this limit) via the continuous 2nd order phase 
transition. 

The corresponding calculations are straightforward for 
the reflecting boundary condition (6). Somewhat unex- 
pectedly from the mathematical standpoint there is only 
one dominant mode (0, 1), which corresponds to the order 
parameter profile $(r) = Acxp(ig*z). All other modes, 
even those unstable being single, are suppressed in the 
presence of this dominant mode. But from the physical 
viewpoint it is quite natural. Indeed, the boundary condi- 
tion (6) does not change the ordering conditions as com- 
pared to that in bulk, where the lamellar morphology is 
the most stable one. Thus, any other ordered morphology 
would cause only some additional frustration without any 
compensation. 

6. The situation changes drastically if there is a prefer- 
ential adsorption of the system particles to the boundary, 



which is naturally described by an additional linear surface 
term [30]: 

F = F bulk + F SUI{ , F surf = (hT/v) J $ (r) dS, (26) 

h being the strength of the preferential adsorption. Obvi- 
ously, the term (26) affects the modes (0, n) with n high 
enough to be of the type (11) and thus favors the coax- 
ial morphologies. Say, the mode (0,2) is dominant in the 
narrow capillaries satisfying the inequality 

P < s , 2 = 3.832 (27) 

if |t| is not too big. If r — > — oo the mode (0, 1) becomes 
dominant since it has the smallest quantitative value of the 
reduced forth vertex. To find the 1st order transition line 
(18) between the modes (0, 1) and (0,2), it is convenient 
to rewrite eq (25) in the reduced variables A = H X / 3 B, 
t = H 2 / 3 t, H 4 / 3 f = 3minF(Si,B 2 ): 

F(B U B 2 ) =bB + t 1 B 1 2 +t 2 B 2 2 (28) 
+ (3/A)Bf + 3B 2 B 2 + (c 22 /A)Bi 

The fields Bi,B 2 correspond to the modes (0,1), (0,2), 
b = 6J (s , 2 )/(A'o, 2 ) 1 / 2 , t 2 = T + Kva/h 2/3 and H = h/R. 

In fact, however, the transition found by minimization 
of the function (28) turns our to be only metastable. The 
case is that, contrary to the slab mode (0, 1), the coaxial 
mode (0, 2) could be compatible (form a mixed state) with 
some helical or angular modes, in particular, simple helix 
(1,1), double helix (2,1), (3,1), (3, 2), (4,1), (4, 2) and so 
on. Thus, to build the phase diagram we are to analyze, 
which of the pure or mixed modes is the most thermody- 
namically favorable and in which temperature interval. In 
fact, it is sufficient to take into account only the modes 
{(1, 1), (2, 1)}, which would start to grow early than the 
coaxial mode (0, 2) in the absence of any preferential ad- 
sorption. The resulting phase diagrams within the inter- 
val (27) are shown in fig. 2. Remember that in the regions 
2b, 3b the long segments with the left-hand helical waves 
alternate with the right-hand ones since the free energies 
of both are identical. 

It is easy to see that increase of the surface field h results 
in increase of the temperature width At (At ~ 10 /i 2 / 3 ) 
of the region, where the helical ordering can be observed. 
Besides, the regions of purely angular ordering 2a and 3a, 
which are just negligible for small h, grow considerably 
with increase of h. It is convenient to visualize the char- 
acter of the appearing angular and/or helical ordering by 
the curves r(4>) satisfying equation 

$(r,0)=O, (29) 

where the order parameter $(r) is given by expression 
(24) with mi = 0,m = 2 and m 2 = 1(2), n 2 = 1. The 
curves (29), which could be considered as intersections of 
the conditional interfaces between the domains with the 
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Fig. 2: The reduced phase diagrams of ordering in the most narrow capillaries satisfying condition (27) for h = 100 (a), h = 1 
(b) and h — 0.01 (c) as well as (d) the dependences of the pitches d (normalized on the lamellar period in bulk D) of the 
ordinary (1) and double (2) helices on the reduced capillary radius p. The regions 1 and 4 correspond to the coaxial (0,2) and 
slab (0, 1) modes, the regions 2 and 3 do to the mixed (0, 2) + (1, 1) and (0, 2) + (2, 1) states, respectively. The dotted horizontal 
and vertical lines on figures (a)-(c) correspond to the bulk ODT temperature r = and the boundary (27), respectively. The 
dotted vertical lines 1' and 2' on figure (d) are the asymptotes p — si t i and p — S2,i of the curves 1 and 2, they correspond to 
the dashed vertical lines on figures (a)-(c) separating the regions a and b of the axially non- modulated and modulated states 
of the same modes. The dotted horizontal line d = D on figure (d) is the common asymptote of the curves 1 and 2 at p — > oo. 



A 




B 



Fig. 3: The curves (29) (domain interfaces) in the 1-thread 
helix (0, 2) + (l, 1) and double DNA-like helix (0, 2) + (2, 1). The 
outer bold circle is the capillary cross-section, the inner one is 
the interface for the pure (0, 2) mode (x — oo). With decrease 
of x the 1-thread helix interface (a) becomes more and more 
eccentric (the plotted dashed curves correspond consecutively 
to x = 5, 3, 2.25, 1.75, 1.5) until the interface (bold solid curve) 
touches the capillary wall at x = 1.445. Next the contact 
area extends (see the thin solid curves corresponding to x = 
1.4, 1.2, 0.8, 0.5, 0.25, 0.1) and at x — the interface is just the 
diameter AB. The double helix interface (6) undergoes a similar 
evolution but the fact that it is centrosymmetrical. 



cross-sections normal to the capillary axes, are plotted in 
fig. 3 for various values of the ratio x = Ai/A 2 , where Ai 
are the coefficients appearing in the expression (24) . 

Beyond the interval (27) the mode (0, 2) is axially mod- 
ulated and does not contribute into the surface term (26) 



anymore. If |t| is not too big, then the modes (0,n) 
are dominant within the intervals so. n -i < P < so.n)- 
With the temperature decrease the pure coaxial modes 
are subsequently replaced by their various combination 
with the compatible helical or angular modes until the 
slab mode (0,1) wins finally for r — > — oo. Thereby, the 
sets of relevant compatible modes are fast growing with in- 
crease of the reduced capillary radius p and, thus, n, which 
makes calculation of the dominance maps rather cumber- 
some. E.g., the relevant sets are {(1, 1), (2, 1), (3, 1), (4, 1)} 
for n = 3 and {(1, 1), (2, 1), (3, 1), (4, 1), (3, 2), (4, 2)} for 
n = 4 (other modes with s m n < So, 3(4) are incompati- 
ble with the corresponding coaxial modes). The occurring 
morphologies appear as ordinary, double (DNA-like) or 
multi-helix structures (sometimes coaxial ones with the 
same pitch) entwining central straight rod. These mor- 
phologies resemble those seen under simulation of the 
asymmetric (with non-zero cubic term) block copolymers 
in refs [1,9] and are rather different from the lamellar mor- 
phology to be observed in the bulk for the case considered 
here (a = 0). 

Summarizing, we have shown that emergence of var- 
ious helical and other complex morphologies in systems 
undergoing ordering when they arc confined to cylindrical 
geometry is a rather general result of interplay between 
the geometrical and energetic effects of confinement even 
in the seemingly trivial case of the reflecting boundary 
condition (6) and zero cubic term. For this case a de- 
tailed "dominance map" (in the limit Gi — > it is just 
the phase diagram) including various complex structures, 
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is presented. 

The actual phase diagram of particular real systems 
could differ from that build in fig. 2. First, the cubic 
term, if present, would contribute to stability of the found 
morphologies formed by straight and curved rod-like units. 
Next, the numerical values of the effective forth vertices 
Cij, which appear in the 2-amplitude Landau expansion 
(25) and determine eventually the shape of the phase 
diagrams shown in fig. 2, depend on the assumption of 
locality of the forth vertex in the Landau Hamiltonian 
(1) used in this paper. If the vertex is non-local and 

4 

reads J F (ri, r 2 , r 3 , r 4 ) J| $ (rj) dr i} then an "angle de- 

i=l 

pendence" [14,31,32] of the forth vertex appears. When 
strong enough, the dependence is known to change the 
sequence of the stable morphologies even in bulk [31-33]. 
Thus, it could change the phase diagram shown in fig. 2 
either. 

Our results suggest also that it could be fruitful to an- 
alyze experimental or numerical data in terms of Bessel 
functions or of their underlying symmetry. More detailed 
presentation of the Landau theory including discussion of 
other possible boundary conditions and microscopic the- 
ory of block copolymer will be given in an extended version 
of this paper. 

* * * 
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REFERENCES 

Wu Y., Cheng C, Katsov K, Sides S.W., Wang 
J., Tang J., Fredrickson G.H., Moscovits M. and 
Stucky CD., Nature Materials, 3 (2004) 816. 
Shin K., Xiang H., Moon S.I., Kim T., McCarthy T.J. 
and Russell T.P., Science, 306 (2004) 76. 
Xiang H., Shin K., Kim T., Moon S.I., McCarthy T.J. 
and Russell T.P., Macro-molecules, 37 (2004) 5660. 
Xiang H., Shin K., Kim T., Moon S.I., McCarthy T.J. 

Macromolecules, 38 (2005) 1055. 
Kim T., Moon S.I., McCarthy T.J. 
Journal of Polymer Science: Part B: 
Polymer Physics, 43 (2005) 3377. 

Yiming Sun, Martin Steinhart, Danilo Zschech, 
Rameshwar Adhikari, Goerg H. Michler and Ulrich 
Gosele, Macromol. Rapid Commun., 26 (2005) 369. 
He X., Song M., Liang H. and Pan C, J. Chem. Phys., 
114 (2001) 10510. 

Feng, J. and Ruckenstein, E., Macromolecules, 39 
(2006) 4899. 

Bin Yu, Pingchuan Sun, Tiehong Chen, Qinghua Jin, 
Datong Ding, Baohui Li and An-Chang Shi, Phys. Rev. 
Lett., 96 (2006) 138306. 
[10] Sevink G.J. A., Zvelindovsky A.V., J. G. E. M. 
Fraaije H. and Huinink P., J. Chem. Phys., 115 (2001) 
8226. 



[2] 
[3] 
[4] 
[5] 

[6] 

[7] 
[8] 
[9] 



and Russell T.P., 
Xiang H., Shin K. 
and Russell T.P. 



[11 

[12; 

[13; 

[14 
[15; 

[ie; 
[17; 

[is; 

[19; 

[20; 

[21 

[22; 

[23; 

[24 
[25 

[26; 
[27; 
[28; 
[29; 



[30 
[31 
[32' 
[33' 



and Garbary R.A., Macro- 



Statistical Physics, 
Soc, 68 (1979) 



Li W., Wickham R.A. 
molecules, 39 (2006) 806. 

Landau L.D. and Lifshitz E.M. 
Part 1 Pergamon Press 1980. 

de Gennes P.G., Faraday Disc. Chem. 
96. 

Leibler L., Macromolecules, 13 (1980) 1602. 
Erukhimovich I. Ya., Polymer Sci. U.S.S.R, 24 (1982) 
2223, 2232. 

Landau L. D., Phys. Zs. Sowjet, 11 (1937) 545. 
Kats E.I., Lebedev V.V. and Muratov A.R., Phys. 
Rep., 228 (1993) 1. 

Brazovskii S.A. and Dmitriev S.G., Sov. Phys. JETP, 
42 (1975) 497. 

Wright D.C. and Mermin N.D., Rev. Mod. Phys., 61 

(1989) 385 

Nabutovskii V.M., Nemov N.A. and Peisakhovich 
Yu.C, Sov. Phys. JETP, 52 (1980) 1111. 

Ho ye J. S. and Stell G., J. Phys. Chem, 94 (1990) 7899. 

Borue V.Yu. and Erukhimovich I.Ya., Macro- 
molecules, 21 (1988) 3240. 

Joanny J.-F. and Leibler L., J. Phys. (Paris), 51 

(1990) 545. 

Ciach A. and Babin V., J.Mol.Liq., 112 (2004) 37. 

Bateman H. and Erdelyi A., Higher Transcendental 
Functions, Vol. 2 Mc Graw-Hill Book Company, inc., New 
York Toronto London 1953 

Angerman H.J., Johner A. and Semenov A.N., 
Macromolecules,^ 2006 6210. 

Stepanow S. and Fedorenko A. A., Phys. Rev. E,73 
2006 031801. 

Kerson Huang, Statistical Mechanics John Wiley and 
Sons, Inc., New York -London 1994 

Grosberg A.Yu. and Khokhlov A.R., Statistical 
Physics of Macromolecules American Insitute of Physics, 
NY, 1994 

Fredrickson G.H., Macromolecules, 20 (1987) 2535. 
Erukhimovich I.Ya., JETP Lett, 63 1996 459. 
Erukhimovich I.Ya., Eur. Phys. J. £,18 2005 383. 
Smirnova Yu.G., G. ten Brinke and Erukhimovich 
I.Ya., J. Chem. Phys., 124 2006 054907. 



p-6 



